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ABSTRACT 

We have used the seven year Wilkinson Microwave Anisotropy Probe (WMAP) data in order to up- 
date the measurements of the intensity signal in the G159.6-18.5 region within the Perseus Molecular 
Complex, and to set constraints on the polarization level of the anomalous microwave emission in the 
frequency range where this emission is dominant. At 23, 33 and 41 GHz, we obtain upper limits on the 
fractional linear polarization of 1 .0, 1.8 and 2.7% respectively (with a 95 per cent confidence level). These 
measurements rule out a sig nificant number of model s based on magnetic dipole emission of grains that 
consist of a simple domain dDraine & Lazarian 1 19991) as responsible of the anomalous emission. When 
combining our results with the measurement obtained with the COSMOSOMAS experiment at 1 1 GHz 
( Battistelli et al.l | 2006l). we find consisten cy with the predictions of the electric dipole and resonance re- 
laxation theory (iLazarian & DrainduOOOl) at this frequency range. 



Subject headings: diffuse radiation - radiation mechanisms: 
ISM:individual(G159.6-18.5) - cosmic microwave background 



general - radio continuum:ISM 



1. Introduction 

The dust-correlated microwave (10-60 GHz) emis- 
sion detected by several cosmic microwave back- 
ground (CMB) ex periments performed in th e two last 
decades (COBE dKogut & Hinshawl Il996l) : OVRO 
at 14.5 and 32 GHz dLeitch et al.1 1 19971) : S askatoon 



1998) from very small (N < 10 atoms) rapidly ro- 



dde Oliveira-Costa et aljl997l): 19 GHz (de Oliveira-Costa et al 



tating (~ 1.5 x 10 s ) carbon based molecules 
in the interstellar medium (the so-called "spinning 

dusf) appears to reproduce well the observational 

I — 1 1 1 1 

const r aints (see e.g. [ Watson et al. 2005; Casassus et al. 

20061: llglesias-Grofhl l2006t iDickinson et al.1 120071 
Tibbs etal l2010h . The detailed theoretical description 



1998b: Tenerife dde Oliveira-Costa et alJI 1999112001), 
COSMOSOM AS dWatson et al.l2005l:lHildebrandt et al 
20071) : VSA dScaife et al.1 120071: iTibbs et al.1 l2010l) : 



)f this family of models h as been recently updated in 



Ali-H aimoud et al. I d2009h . also including the effect of 



and references therein) suggests the existence of a 
new continuum microwave emission mechanism un- 
like to the three well-known Galactic mechanisms: 
synchrotron, free-free and thermal dust emission. A 
remarkable observational effort has been devoted to 
the understanding of the intensity and polarization 
properties of this "anomalous" microwave emission, 
among other reasons, because of the importance of an 
accurate foreground correction of the CMB maps at 
low frequencies. 

Among the various scenarios prop osed to explain 
this emission, electric dipole radiation ( Draine & Lazarianl 



rotation of the dust grains around a non-principal axis 
dSilsbee et alJl20ld) . 



An alternative explanation based on magnetic 
dipole emission of spinning dus t has also been pro- 
posed ( Draine & Lazarian 19991) . Measurements of 
the polarization properties of the anomalous mi- 
crowave emission may potentially distinguish between 
these two models. According to ILazarian & Draine 
d2000h . electric dipole radiation from spinning dust 
would be polarized at low frequencies, reaching a 
maximum (6-7%) at 2-3 GHz and dropping to 4-5% 
at 10 GHz and progressively decreasing at higher fre- 
quencies. Polarization from magnetic dipole emission 
predicts a different frequency behaviour and stronger 
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G1 59.6-18.5 



linear polarization depending on the composition and 
shape of the emitting particles. 

From the observational point of view, there is lit- 
tle information in the literature about the polarization 



proper ties of the anomalous emission. iKogut et al 
J2007I) used the full-sky WMAP 3-year data to con- 
strain the polarization fraction of a diffuse anoma- 
lous component traced by the dust morphology. They 
concluded that the polarized anomalous emission con- 
tributes less than 1 % of the observed polarization sig- 
nal variance in any of the five WMAP bands. 

Only few attempts have been made to determine 
the polarization of the anomalous microwave emis- 
sion in individual objects, and in the frequency range 
where this contribution is dominant (10-50 GHz) 



Battistelli etall d2006l) . using data from the COSMO- 



SOMAS experiment on the Perseus molecular com 
plex, reported at 11 GHz n = 3.4 + }-g% (95% confi 



dence level). More recently Mason et al. ( 20091) . using 
the Green Bank Telescope at 9 GHz, obtained an up- 
per limit of linear polarization of 88 fiK at 95.4% of 
confidence in the Lynds 1622 dark cloud. Both the 
Perseus molecular complex and the Lynds 1622 cloud 
are regions where anomalous microwave emission ap- 
pears to dominate on other emission processes in the 
frequency range 10-50 GHz. 

In this work, we present new measurements of the 
polarization of the anomalous microwave emission in 
G159.6-18.5 (within the Perseus molecular complex) 
using the WMAP (Wilkinson Microwave Anisotropy 
Probe) 7-year data. Although the results only provide 
upper limits to the polarized emission in the region, 
they still constitute a strong constraint on the physical 
mechanism responsible for the emission. 

2. The G159.6-18.5 region 

The Perseus molecular complex is a giant molecular 



cloud located at a distance of 260 pc (ICernicharo et al 



19851) . Our region of interest is G159.6-18.5 , a dust 
feature in this molecular complex, observed in the 
IRIS0 maps, which appears as a slightly broken ring 
with a diameter * 1.5°, and a intensity of the order of 
100-200 MJy sr 1 at 100 fim and 5-10 MJy sr" 1 at 
12 jum (see FigureHJ. 

Originally, G15 9.6-18.5 was considered to be a su- 
pernovae remnant ( Pauls & Schwartz! 19891: Fiedler et al 




RIGHT ASCENSION (J2000) 



'IRIS are improved version of IR AS maps, for details see 
iMiville- Deschenes & Lagachd fe005l) . 



Fig. 1. — : IRIS images at 12/im, 60/xm (middle) and 100/im 
(bottom panel) of the G159.6-18.5 region, with a field of 
view of 4.5° x 4.5°, and centered in R.A. (12000)= 55.4° and 
Dec.= +31.8°. All maps are in units of MIy sr -1 . 
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1994), but later observations showed that this is a Hn 
region driven by the O9.5-B0 V st ar HP 278942, 



locate d at the center of the ring dd e Zeeuw et al. 



19991). More rec ent studies dAndersson et all l200oT 



Ridg e etal.|[2006h ) concluded that G159.6-18.5 was 
indeed an expanding H n bubble that has emerged from 
the outer edge of the cloud. 

hereafter W05) carried out 



Watson et all (12005 



a detailed study of the spectral energy distribution 
(SED) of G159.6-18.5 in a wide frequency range, 
from 408 MHz to 3000 GHz. This study combined, 
among others, the observations performed with the 
COSMOSOMAS experiment at 11, 13, 15 and 17 GH z 



Fernandez-Cerezo et al. I l2006l) . 
and the 3rd year WMAP data dBennettetalJl2003l) . 



(Gall egos et all 12001 



The total emitting region in G159.6-18.5 was slightly 
resolved in COSMOSOMAS and WMAP data, and 
was modelled using an ellipse with 1.6° x 1.0° and 
P.A.= 51°, centered on R.A.= 55.4° and Dec.= 
+31.8°. The SED for that region showed clear ev- 
idence for anomalous microwave emission, with a 
peaked spectrum around 20-30 GHz, indicative of 
spinning dust. Indeed, W05 showed that an adequate 
fit to the SED can be achieved from 10 to 50 GHz 
only when including three components in the analysis: 
(a) optically thin free-free emission; (b) vibrational 
dust emission with 7d us t = 19 K and emissivity index 
of 1.55; and (c) a combination of the spinning dust 
models of Draine & Lazariar] ( 1998 ) for warm neutral 
medium and molecular cloud (0.8 WNM + 0.3MC). No 
bright unresolved source that could be ultra-compact 
H n region or gigahertz -peaked source could be found. 

Detailed observations of G159.6-18.5 with the 
Very Small Array (VSA) interferometer at 33 GHz 



•ay 

dTibbs et al.1 1201 On and an angular resolution of 10- 



40 arcmin, showed that the region consists of five dis- 
tinct components, all of which are found to exhibit a 
emission at 33 GHz which is highly correlated with 
the far-infrared emission. The most interesting result 
is that the VSA resolved out most of the emission in 
the region, as those five components contribute to only 
« 10 per cent to the total flux density of the diffuse ex- 
tended emission detected in W05. Therefore, the bulk 
of the anomalous emission in G159.6-18.5 is diffuse. 

Concerning the polarization level of the anomalous 
microwave emission in the region, the only measure- 
ment at the relevant fr equencies ( 1 0-30 GHz) done 
so far was presented in iBattistelli et al.l ( 120061) . Us- 
ing dual orthogonal polarizations with the COSMO- 
SOMAS experiment, the resulting total polarization 



level in G159.6-18.5 at 11 GHz was found to be n = 
3.4^[ -g% (95% confidence level). Based on this value, 
they concluded that this weak detection of polarization 
would be associated to the sp inning dust grains. Re- 
cently, [§iich]fcReich| |200^ suggested that G159.6- 
18.5 is acting as a Faraday screen hosting a strong 
regular magnetic field, which would rotate the back- 
ground polarized emission. They suggest that a Fara- 
day screen model with a rotation measure of RM = 
190 rad m -2 would explain the 11 GHz observations, 
as well as their new 1 1 cm Effelsberg observations in 
the re gion (see Fig. 6 i n that paper). In any case, even 
if the Battistelli et al. ( 2006h result is only an upper 
limit to the polarization level of the anomalous emis- 
sion, its amplitude favours the electric dipole emission 
model with resonance p aramagnetic relaxation (see 
Lazarian & Draine l2000l) ). 



In this context, the seven-year WMAP data rep- 
resent an opportunity to update the measurements of 
intensity signal in the Perseus region, as well as to 
constrain the polarization level of the anomalous mi- 
crowave emission in the frequency range where this 
emission is domi nant (20-30 GHz). As shown by 
iTibbs et alJ (120101) . 90% of the emission in the region 
is diffuse, so the angular resolution of WMAP data is 
sufficient to provide a reliable measurement. 

3. Data and Methodology 

3.1. WMAP data 

In this paper, we use the seven-year WMAP data 
products ( Jarosik et al. 2010l) . which are available in 
the LA MBDA 2 _ website in t he HEALPdJI pixelisation 
scheme (iGorski et al.ll2005l) . In particular, we use the 
smoothed I, Q and U maps, for each of the five fre- 
quency bands centered at 23, 33, 41, 61 and 94 GHz 
(K, Ka, Q, V and W respectively). The original an- 
gular resolution of these maps is approximately 0.82, 
0.62, 0.49, 0.33, and 0.21 deg respectively, but they 
are degraded to a common resolution of 1 deg. For 
each map, a full description of the noise covariances 
(II, QQ, UU, QU) is provided. This information will 
be used in our analyses. We would like to stress that 
the WMAP polarization maps are derived from the dif- 
ference of two differential measurements. Thus, the 
WMAP measures a double difference in polarized in- 
tensity, not the intensity of the difference of the electric 



http ://lambda.gsf c . nasa. gov/ 
3 http://healpix.jpl. nasa.gov 
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Fig. 2. — : Neighborhood of the Perseus Molecular Complex as seen by WMAP-7 at 33 GHz (Ka band), with field of view of 
23° x 23° centered in the G159.6-18.5 region. All images (I, right panel; Q, center; and U, left panel) are smoothed to 1° resolution, 
and in units of mK (thermodynamic). The black circles define the circular apertures around the objects studied in this paper (see 
text for details). 



field as with interferometers and correlation receivers 



(for further details. seelPage et al.l 120071 : iJarosik et al 
2007l:lKogutetal.ll2003f) . 



From those Q and U maps at each frequency band, 
we can obtain the maps of the polarization intensity as 
P = V2 2 + while the angle between the polar- 
ization direction of the electric field and the Galactic 
meridian can be obtained as y — I arctan({//0. 

Figure [2] shows a detailed view of the Perseus 
Molecular Complex as seen in the Ka-band of WMAP 
(33 GHz). In this map, two strong sources are clearly 
visible: G159.6-18.5 (placed at the center of the im- 
age), and the Hn region NGC1499. This later ob- 
ject will be used as null hypothesis in our study, as 
free-free emission is known to be un polarized. Th is 



null-te st was also used in the study of iBattistelli et al 



(120061) . Finally, the figure also displays six nearby re- 
gions which have been selected due to their low emis- 
sion in intensity. These regions will be used below to 
characterize the behaviour of the background diffuse 
emission in the surroundings of G159.6-18.5 . 

3.2. Ancillary data 

3.2.1. Intensity 

As a reference for the SED in G159. 6-18.5 , in 
this work we use the data points from Table 1 in 
W05. In that table, in addition to the COSMOSOMAS 
and WMAP me asurements, there are also data points 
at 0.408 MHz dHaslam et all Il982h : at 1420 MHz 



dReich&Reichlll988l) : and at 100, 140 and 240 /an, 
based on the DIRBE maps dSilverberg et alj|l993l) . 

Here we have re-evaluated the intensity measure- 
ments on the Haslam, WMAP and DIRBE maps, using 
the methodology described in Section [3~3l The derived 
fluxes are slightly higher than those presented in W05 
(by a factor of 5-13 per cent), although as pointed out 
below, this is probably due to the different methodol- 
ogy used for the flux determination in W05. 

After W05. ISemenova et alJ d2009t) used the RATAN- 
600 telescope to obtain independent measurements in 
the region. Although an estimate of the error bar of 
these new measurements is not provided in their pa- 
per, the overall shape of the SED in G159.6-18.5 in the 
frequency range 1-20 GHz is fully compatible with the 
results of W05, showing a rising spectrum in that fre- 
quency range. 

Recently, Rei ch & Reich d2009h presented prelimi- 
nary results of the emission in G159.6-18.5 at 11 cm 
and 6 cm with the Effelsberg and the Urumqi tele- 
scopes, respectively. These results support that the 
emission in the frequency range 2.7-5 GHz is compati- 
ble with an optically thin thermal gas, as expected from 
theSEDinW05. 

3.2.2. Polarization 

The only reported measurement to date of the ra- 
dio polarized intensity in G159.6-18.5 c omes from the 
COSM OSOMAS experiment at 1 1 GHz dBattistelli et all 
20061) . The measured fractional polarization is n = 
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Fig. 3. — : Close view of G159.6-18.5 as seen by WMAP-7 at 33 GHz (Ka band), with a field of view of approximately 8° x 8°. The 
images show I (right frame), Q (central frame) and U (left frame) at 1° resolution, and in units of mK (thermodynamic). The small 
aperture at the northwest of the center of G159.6-18.5 corresponds to the solid angle defining the quasar 4C+32.14. The pixels 
within this small aperture were not considered in our flux estimation. 



3A + _[ 5 g % (95% C.L.). 



Recently, iReich & Reich (2009) presented prelim- 
inary results of the polarized emission in G159.6- 
18.5 at 11 cm (~ 2.7 GHz) with the Effelsberg tele- 
scope. They reported a significant degree of polariza- 
tion along the ring, which they interpreted as an indica- 
tion that G159.6-18.5 acts as a Faraday screen rotating 
the polarization angle of the background emission in 
the region. However, this effect will be strongly sup- 
pressed at the peak frequencies of the anomalous emis- 
sion (20-30 GHz) as discussed below. 



Kogut et alJ (120071) presented a full-sky model of 



the polarized Galactic microwave emission, consider- 
ing in their analysis the polarization from the thermal 
dust and synchrotron emissions. The resulting polar- 
ized dust model obtained in this paper shows that the 
Perseus molecular complex lies near to the boundaries 



of th e P06 mask (see mask definition in iPage et al 



20071) . where the expected dust fractional polariza- 
tion (U dust ) reaches 6.4 + 1.3% and 2.0 + 1.4% for 
\b\ < lOdeg and \b\ > lOdeg, respectively, inside 
P06; and 3.6 ± 1.1% outside the P06 mask. This 
fractional degree of polarization will be of importance 
when studying the high frequency behaviour of the po- 
larized emission in the G159.6-18.5 region. 

3.3. Determination of the Spectral Energy Distri- 
bution in intensity and polarization 

Figures[2]and[3]show that there is no obvious struc- 
ture seen within the G159.6-18.5 region in the polar- 



ization maps. In this case, we can not apply a similar 
method to the one used in W05 for the flux determina- 
tion. In that paper, a simple elliptical Gaussian model 
was used to describe the region, and fluxes were ob- 
tained using a direct fitting to the data. 

Here, we will adopt a direct aperture integration, 
or "ring analysis", both for the intensity and po- 
larization determinations. The practical implemen- 
tation of thi s method has been previously used in 
other works dBennett et alJll993b iBandav et al1|l996t 



Hernandez-Monteagudo & Rubino-Martrnl2004 ). The 
basic idea is as follows. For a certain frequency map, 
we define a circular aperture around the object with ra- 
dius r\ , defining a solid angle f2s around the source. In 
order to correct for the background emission, a circular 
corona with radius r\ < r < rj_ is used. In order to pre- 
serve a similar area both for the source and the back- 
ground flux estimation, we use the relation ri — r\ a/2 
between the two apertures. This method provides an 
efficient correction of the background emission around 
the source for apertures (ri values) comparable to the 
beam size, although for larger apertures might be bi- 
ased due to the background or CMB contributions. 
The error bar associated to this measurement contains 
the quadratic sum of two contributions, the pure in- 
strumental noise and the increase of the variance due 
to the CMB fluctuations. It is important to note that 
the evaluation of both contributions should account 
for correlation terms. In the case of the CMB, this is 
introduced by means of the two-point correlation func- 
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tion, which is evaluated using the (measured) WMAP 
angular power spectrum. In the case of the noise, the 
correlation between pixels appears after the smoothing 
process, and is taken into account in a pixel-by-pixel 
basis. 

As we will see below, in the case of the selected 
aperture for G159.6-18.5 (r\ = 2°) and the instrumen- 
tal noise level of the WMAP data, the CMB contri- 
bution to the error bar is dominant for the intensity 
Stokes parameter, while its contribution is negligible 
for Stokes Q and U. Although the main focus of this 
paper is the polarization determination, for illustration 
we have also considered how the intensity measure- 
ments change with our aperture method when a CMB 
correction is performed before the flux determination. 
For that discussion, we have used the Internal Lin- 
ear Combination (ILC) map pr ovided by the WM AP 
team in the LAMBDA webpage dGold et al.l201dh . Al- 
though this map does not provide an optimal separa- 
tion of CMB component, it provides a reasonable ap- 
proximation which can be used to estimate its impact 
on the final intensity measurements. 

4. Results and Discussion 
4.1. Control regions 

4.1.1. Characterization of the background emission 

G159.6-18.5 is relatively close to the galactic plane, 
in a region with a non-negligible diffuse emission 
which varies markedly with the position. In order to 
characterize the behaviour of this background emis- 
sion, we have selected six regions with low emission 
in intensity in the surroundings of the Perseus molecu- 
lar complex. The central coordinates of the six regions 
are shown in table Q] Around each of them, we have 
chosen the same aperture as for G159.6-18.5 (n = 2°) 
for the flux determination. 

Table|2]summarizes our results. In addition to the I, 
Q, and U determinations at each frequency, we also in- 
clude in the last two columns the values for the fluxes 

Table 1:: Identification number and coordinates of the six 
control regions in the surroundings of G159.6-18.5 used for 
the background characterization. The aperture used in all 
cases and at all frequencies is r\ =2°. 



ID 


Al 


A2 


A3 


A4 


A5 


A6 


r.a. n 

DEC. [°] 


52.19 
36.88 


47.30 
35.74 


50.21 
28.56 


46.88 
24.34 


60.39 
23.88 


59.70 
27.48 



inside the inner aperture (ri) with no background cor- 
rection. These numbers, quoted as Qi, ac k and Uback, 
will be used below for the analysis of the Faraday 
rotation in the region (Sect. 14.2b . The table also in- 
cludes an estimation of the total linear polarization 
P = -\/Q 2 + U 2 , and also the upper limit on this lin- 
ear polarization Po that we would obtain in the null- 
hypothesis case using a zero signal (Q = U = 0) and 
the same noise levels. As in all these values no detec- 
tion is found, we quote upper limits which correspond 
to the 95% confidence level derived from a maximum 
likelihood approach based on the Q and U measure- 
ments. 

We first note that at all frequencies, the main con- 
tribution to the error budget in intensity comes from 
the CMB part, which for illustration purposes is writ- 
ten separately in the table [2] (values in brackets in the 
second column). As one would expect, its relative con- 
tribution rises with approximately v 2 dependence, as 
the temperature sensitivity of WMAP is comparable in 
all channels. In the case of Q and U, the CMB contri- 
bution is negligible and it is not included in the error 
bar. 

Finally, as a summary of the results in table [2] we 
present in table [3] the peak-to-peak variations (l p - p ), 
means and the root-mean-squares of these values for 
the six regions. We note that in the case of intensity, 
the variations in the three first frequencies (23, 32 and 
41 GHz) is even larger than the expected CMB contri- 
bution, which implies that in addition there is a signif- 
icant contribution to the error budget coming from the 
variations of the diffuse background in the surround- 
ings. This effect is not seen in polarization, so in this 
case the error bars seem to be dominated by the in- 
strumental noise contribution. However, we note that 
the average values of Qback and Uback show a small de- 
gree of diffuse polarization in the background, which 
is very well corrected by the ring analysis. 

4.1.2. NGC 1499: Null-hypothesis in polarization 

The California Nebula is a Hn region close to 
G159.6-18.5 (its optical coordinates are RA(J2000)= 
04 /, 03" , 18" andDec= +36°25.3'). Since the dominant 
emission mechanism in diffuse H n regions at these fre- 
quencies is the free-free emission, they are expected to 
be practically unpolarized, so NGC 1499 will be our 
null-hypothesis for the polarization. 

California is also an extended object at this res- 
olution. For consistency, we use the same aperture 
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Table 2:: Results obtained for the six control regions in the surroundings of G159.6-18.5. Stokes I, Q, U are measured using 
the "ring analysis" method with an aperture of rj = 2° at all frequencies. Error bars in all cases include the instrumental noise 
contribution. In the case of Stokes I, the number in brackets indicates the expected CMB contribution to the error budget, which 
has to be added in quadrature to the instrumental noise. Columns 5 and 6 present the derived upper limits (at the 95% confidence 
level ) on the linear polarization P, as well as upper limit P that would be obtained in the null-case (i.e. using Q=U= 0). Last two 
columns (Qtack and LWa) show the direct flux determination within the aperture r\ without the diffuse background correction (i.e. 
without the subtraction of the flux in the circular corona around n). 



Region Al 



y 


[Lrrlzj 




i UyJ 




T T rTi*l 


" UyJ 


ro UyJ 


\lback [Jyl 


Uftor* UyJ 




23 


-3.57 


± 0.13 (± 1.49) 


-0.27± 0.16 


-0.27± 0.21 


< 0.72 


< 0.46 


1.12± 0.13 


-1.28+ 0.16 




33 


-6.05 


± 0.23 (± 3.00) 


-0.39± 0.30 


0.43± 0.38 


< 1.48 


< 0.82 


0.39± 0.26 


-0.85± 0.33 




41 


-7.71 


±0.32 (±4.59) 


-0.02± 0.41 


-0.13±0.54 


< 2.10 


< 1.17 


0.02± 0.37 


-0.86± 0.49 




61 


-13.86 


± 0.80 (± 9.67) 


0.09± 1.02 


0.16± 1.31 


< 3.29 


< 2.86 


0.93± 0.95 


-0.97± 1.22 




94 


-20.81 


± 1.92 (± 19.98) 


-3.98±2.44 


0.95±3.19 


<7.26 


<6.95 


-1.00± 2.29 


9.74± 3.00 


Re; 


jion A2 




















23 


3.97 


± 0.13 (± 1.50) 


0.23±0.17 


0.26± 0.21 


<0.70 


<0.46 


1.75+ 0.13 


-1.30+ 0.17 




33 


5.91 


± 0.23 (± 3.02) 


0.01±0.30 


0.04± 0.38 


< 1.09 


<0.84 


0.68± 0.26 


-1.11+ 0.33 




41 


7.43 


±0.34 (±4.61) 


-0.26± 0.44 


0.09± 0.56 


< 1.37 


< 1.23 


0.36± 0.40 


-1.33+ 0.51 




61 


13.58 


±0.81 (± 9.72) 


-0.78± 1.03 


-0.81+ 1.35 


< 3.33 


<2.93 


0.78± 0.95 


-0.10± 1.25 




94 


28.66 


± 1.98 (± 20.09) 


0.40± 2.50 


-4.84± 3.29 


<7.92 


< 7.16 


5.48± 2.34 


2.54±3.10 


Re; 


jion A3 




















23 


6.14 


±0.13(± 1.49) 


-0.25±0.17 


-0.60± 0.23 


< 1.02 


<0.49 


1.03+ 0.13 


-1.74± 0.18 




33 


5.59 


± 0.24 (± 2.99) 


-0.20±0.31 


0.47±0.41 


< 1.35 


<0.89 


0.32± 0.27 


-0.52± 0.35 




41 


6.40 


± 0.36 (± 4.57) 


-0.06± 0.46 


-0.07± 0.59 


< 1.79 


< 1.29 


0.42± 0.41 


-1.07± 0.54 




61 


8.86 


± 0.89 (± 9.63) 


0.48± 1.13 


1.22+ 1.48 


<4.54 


< 3.21 


1.3 1+ 1.05 


0.52± 1.39 




94 


14.65 


±2.17(± 19.91) 


0.41±2.75 


0.22± 3.59 


<9.90 


< 7.81 


-1.84+ 2.59 


5.91+ 3.39 


Re: 


gion A4 




















23 


-2.90 


±0.14(± 1.49) 


0.00±0.18 


-0.74± 0.24 


< 1.14 


< 0.51 


1.14+ 0.14 


-1.51+ 0.19 




33 


-3.41 


± 0.26 (± 2.99) 


0.18±0.33 


-0.36± 0.43 


< 1.54 


<0.93 


0.72± 0.28 


-0.17+ 0.37 




41 


-4.45 


±0.41 (±4.58) 


0.74± 0.50 


1.35±0.75 


<2.40 


< 1.58 


1.09± 0.46 


0.79± 0.68 




61 


-7.20 


±0.91 (± 9.65) 


-0.03± 1.15 


-2.88± 1.51 


<4.79 


< 3.28 


0.93± 1.07 


-2.89± 1.40 




94 


-11.37 


±2.26 (± 19.94) 


1.06± 2.88 


1.03+ 3.71 


<9.35 


< 8.12 


5.22± 2.71 


12.62± 3.48 


Region A5 




23 


-2.02 


±0.13(± 1.50) 


0.22±0.17 


-0.21±0.22 


<0.66 


<0.48 


1.39+ 0.14 


0.09± 0.17 




33 


-1.65 


± 0.24 (± 3.02) 


0.25±0.31 


-0.47± 0.40 


<0.96 


<0.87 


0.73± 0.27 


-0.29± 0.35 




41 


-1.83 


±0.36 (±4.61) 


0.90± 0.48 


-0.26± 0.57 


<2.04 


< 1.28 


1.07± 0.43 


-0.1 0± 0.51 




61 


-2.10 


± 0.86 (± 9.72) 


-0.31+ 1.11 


0.20± 1.38 


<4.95 


< 3.05 


-0.02± 1.03 


1.20± 1.28 




94 


-11.06 


±2.12(± 20.08) 


1.81±2.74 


-2.49± 3.38 


< 8.07 


<7.50 


7.06± 2.59 


2.43±3.19 


Region A6 




23 


-3.95 


±0.13(± 1.49) 


0.15+ 0.17 


-0.13+ 0.22 


<0.57 


<0.48 


1.46± 0.13 


-0.47± 0.17 




33 


-2.91 


± 0.24 (± 3.00) 


-0.05± 0.30 


-0.28± 0.40 


< 1.17 


<0.86 


0.80± 0.26 


-0.39± 0.34 




41 


-2.48 


± 0.35 (± 4.59) 


-0.42± 0.45 


0.94± 0.56 


< 1.56 


< 1.25 


0.65± 0.40 


1.16+ 0.51 




61 


-1.34 


±0.83 (± 9.68) 


0.17± 1.06 


-0.09± 1.37 


< 3.44 


<2.99 


-0.26± 0.98 


1.31+ 1.28 




94 


-5.80 


± 2.04 (± 20.00) 


1.73+ 2.59 


1.37+ 3.36 


<7.84 


<7.36 


4.34± 2.44 


5.86+ 3.17 
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Table 3:: Summary of the statistics for the six control regions in table [2] We present the peak-to-peak variations (I p - P ), means 
and standard deviations of the Stokes I, Q and U, together with the mean of the Q/, fl „t and values for each of the five WMAP 
frequencies. All quoted values are given in Jy. 



v [GHz] 




<D 


(I > 


Qp-p 


<Q> 


<Q > ' 


V P - P 


<u> 


<uV ; - 


(Qback) [Jy] 


OW>[Jy] 


23 


10.10 


-0.39 


4.32 


0.50 


0.01 


0.23 


1.00 


-0.28 


0.36 


1.32 


-1.04 


32 


11.96 


-0.42 


4.99 


0.64 


-0.03 


0.24 


0.95 


-0.03 


0.41 


0.61 


-0.55 


41 


15.13 


-0.44 


6.06 


1.32 


0.15 


0.54 


1.61 


0.32 


0.66 


0.60 


-0.24 


61 


27.44 


-0.34 


10.12 


1.26 


-0.06 


0.44 


4.10 


-0.37 


1.39 


0.61 


-0.15 


93 


49.46 


-0.96 


18.70 


5.80 


0.24 


2.16 


6.21 


-0.63 


2.50 


3.21 


6.52 



r\ = 2° for all frequency bands, which after visual 
inspection, is found to fit well with the solid angle 
that defines the source. The central coordinates of the 
emission used for the flux integration are RA= 60.40° 
and Dec= +36.3675°, practically coinciding with the 
peak emission in the 33 GHz map. Our results on the 
I, Q, U measurements, and the derived upper limits on 
the linear polarizatiorQ P and II, as well as the null- 
case Pq and IIo, are presented in table @] 

A simple extrapolation of the DIRBE fluxes to the 
WMAP frequency range shows that the thermal dust 
emission is expected to be negligible between 20- 
100 GHz. Thus, we expect that the main emission 
mechanism is the free-free emission. To check this 
prediction, we have derived the spectral index that 
best-fits the data obtained in the range 23-41 GHz, us- 
ing a power law model (I oc y 8 ). We find /? = 0.08 ± 
0.10, which is apparently in contradiction with the ex- 
pected frequency behaviour for an optically thin free- 
free emitting region at these frequencies (ft as -0.1). 
However, this seems to be a effect of the residual CMB 
contribution to these measurements. If we correct the 
WMAP maps from the CMB emission by subtracting 
the ILC map, the flux estimates at 23, 32 and 41 GHz 
are now found to be 61.6, 59.6 and 59.5 Jy, respec- 
tively. Using these values, we find /3 = -0.07 + 0.11, 
which is now compatible with the expected value. 

Concerning the polarization measurements, at all 
frequencies the Stokes Q and U parameters are found 
to be compatible with a zero level. The upper lim- 
its on P, and n at the 95% confidence level are simi- 
lar to those obtained for the null-case hypothesis (Pq 
and no). Indeed, the full posterior distributions for 
these parameters are found to be very similar, imply- 
ing a non-detection of polarization in the region, as one 
would expect for an H n region. 



4 IT is the fractional linear polarization, denned as If = 100 P/I. 



4.2. Perseus 

Figure [3] shows a detailed view of the G159.6-18.5 
region. The circular aperture used for the flux deter- 
mination corresponds to a radius of r\ = 2°, and is 
centered in RA(J2000)= 55.4° and Dec= 31.8°. Since 
G159.6-18.5 is a extended object, we have used the 
same aperture r\ at all frequency bands. The sub- 
set of pixels within the small aperture at the north- 
west of the center of G159.6-18.5 were removed for 
the flux determination, as they correspond to the loca- 
tion of the quasar 4C+32.14. This small aperture uses 
7"4c+32.i4 = 0.25° at all frequencies, and it is centered 
at RA(J2000)= 54.125° and Dec= 32.308°. Table 
summarizes the results obtained for the region. 

4.2.1. Intensity Signal 

The second column in Table [5] presents the results 
for the intensity measurements in the G159. 6-18.5 re- 
gion. A comparison with the results in W05 shows 
that our fluxes are between 5 and 13 per cent higher 
in the first four channels (23 to 61 GHz), and 10 per 
cent smaller in the last channel (94 GHz). This small 
discrepancy is probably due to the different method- 
ologies used to derive the fluxes. In W05, an elliptical 
model for the region was adopted, so any departure of 
the actual shape of the region from this model might 
introduce these changes. 

Figure|4]shows our results for the intensity signal in 
WMAP and DIRBE maps, together with the ancillary 
data described in section [3~2l Following W05, using 
the W band and the DIRBE data we can fit the ther- 
mal dust emission with a modified blackbody function 
hust K v^"" +2 B(v, Tdus,), where Ta ust is the dust tem- 
perature and Pdust the dust emissivity index. The black 
dashed line represents the best-fit model of the thermal 
dust emission, which has parameters fi^ust = 1 -55 and 
Tdust = 19 K, in full agreement with W05. These val- 
ues are expected for dust in a 'warm neutral medium' 
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Table 4:: Summary of flux measurements for NGC1499. We use an aperture of r\ = 2° at all frequencies to integrate the fluxes for 
I, Q and U Stokes parameters. The upper limits on P and IT are derived using a maximum likelihood analysis. Similarly, we get the 
upper limits Po and TIo for the null-case hypothesis (Q=U=0). As in Table|2] for the intensity column the error budget is separated 
in two contributions, the instrumental noise and the CMB error (within parenthesis). 



v [GHz] 




i[Jy] 


Q[Jy] 


U[Jy] 


P[Jy] 


n [%] 


Po [Jy] 


n [%] 


23 


64.65 


±0.13(± 1.50) 


0.28±0.16 


-0.12+0.21 


<0.63 


<0.98 


<0.46 


<0.71 


33 


65.70 


± 0.23(+ 3.02) 


-0.22± 0.30 


0.44± 0.37 


<0.90 


< 1.40 


<0.82 


< 1.27 


41 


68.92 


± 0.32(± 4.62) 


0.36+ 0.41 


0.15+0.54 


< 1.18 


< 1.84 


< 1.17 


< 1.84 


61 


70.73 


±0.79(+9.73) 


0.58+ 1.01 


0.21± 1.31 


<3.05 


<4.73 


<2.84 


<4.43 


94 


65.61 


± 1.93(± 20.11) 


2.64+ 2.45 


-1.55+3.19 


<7.10 


< 10.90 


<6.91 


< 10.70 



Table 5:: Summary of flux measurements for the G159.6-18.5 region. The meaning of the different columns is the same as in 
TableH 



v [GHz] 




I [Jy] 




Q[Jy] 


U[Jy] 


P[Jy] 


n [%] 


Po [Jy] 


n [%] 


23 


47.76 


±0.13(± 


1.44) 


-0.07+0.16 


0.01+0.22 


<0.48 


< 1.01 


<0.47 


<0.98 


33 


42.97 


± 0.23(± 


2.90) 


0.05+ 0.29 


0.43+ 0.39 


<0.86 


< 1.79 


<0.85 


< 1.76 


41 


37.88 


±0.33(± 


4.44) 


-0.09+ 0.42 


-0.41+0.55 


< 1.29 


<2.69 


< 1.19 


<2.48 


61 


36.36 


± 0.80(± 


9.35) 


0.06+ 1.01 


0.67± 1.33 


<3.44 


<7.23 


<2.88 


<6.03 


94 


70.18 


+ 1.94(+ 


19.32) 


1.13+2.46 


-2.06+ 3.24 


<7.52 


< 15.64 


<7.04 


< 14.65 



WNM). 

The low frequency behaviour is fitted by the dotted 
line to a free-free emission, using the measurement at 
1420 MHz (7.3 + 2.0 Jy) and the typical free-free spec- 
tral index /3yy = -0.12 for an optically thin region. 

As in W05, the extrapolations of the thermal dust 
and free-free emissions can not explain the rising spec- 
trum between 10 and 23 GHz, and a third compo- 
nent has to bee included in the fit. Following W05, 
we co nsider a model based on a linear combination 
of the Draine & Lazariar] ( 1998 ) models for Warm 
Neutral Medium (WNM) and Molecular Cloud (MC). 
Using T3000 = 8.45 x 10 -4 (obtained from the fit 
thermal dust emission) and the ca nonical factor of 
2.13 x 10 24 Hem" 2 = Itiqo (from iFinkbeiner etaT 



20041) . we obtain the column number density of hydro- 
gen atoms N{H) = 1.8 x 10 21 cnT 2 , which differs from 
the value obtained by W05 (1.3 x 10 22 cm~ 2 ) due to 
the dependence of N(H) on the solid angle. Finally, 
considering our N(H) value and using the solid an- 
gle Q. s that defines the source, we find that the best-fit 
to the observed data is obtained with the combination 
l a ~ 0.85WNM + 0.35MC (black dot-dashed line in 
Figure 0). The red solid line in Fig. |4] corresponds to 
the co-added spectral energy distribution model \ tot for 
the three components (\, , = l a + \d ust +!//)• 



4.2.2. Polarization Signal 

Table summarizes our polarization measure- 
ments, which show that the polarized intensity in 
G159.6-18.5 is compatible with zero at the five fre- 
quency bands. Because of this reason, all quoted val- 
ues in columns 5 to 8 correspond to upper limits at the 
95 per cent confidence level. Nevertheless, these upper 
limits on the linear polarization fraction still provide 
extremely important constraints on the nature of the 
anomalous emission. 

Figure [5] represents our upper limi ts from WMAP 
data, t ogether with the 1 1 GHz result of lBattistelli et al 
(12006b . For comparison, we also include the expected 
dust fractional polarization for three cases, namely 
Rdust equal to 6.0% (solid black line), 3.6% (solid 
blue line) and 2.0% (solid red line). This range of 
values for the polarized dust emission in t he region are 



taken t o be consistent with the estimates in lKogut et al 



(120071) . The conclusion is that our upper limits on the 
fractional linear polarization in G159.6-18.5 at these 
microwave wavelengths can be directly translated into 
upper limits on the anomalous emission, at least in 
the frequency range 23-41 GHz, while the dust con- 
tribution has to be taken into account at higher fre- 
quencies. Thus, the upper limits at 23 and 32 GHz 
of n(23 GHz) < 1.0% and n(32 GHz) < 1.8% (95% 
C.L.) constitute the most stringent constraints on the 
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polarization of this emission. 
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Fig. 4. — : Spectral Energy Distribution (SED) for the 
G159.6-18.5 region. The blue squares and black diamond 
correspond to measurements obtained in this paper (WMAP 
and DIRBE bands) and the ancillary data used (COSMOSO- 
MAS and HASLAM bands), respectively. The SED is fit- 
ted with three contributions: i) the free-free emission (I//), 
modelled as a power law with spectral index /}// = -0.12; 
ii) the thermal dust emission, modelled as a modified black 
body (Jdust) with pj lts , = 1.55 and T dust = 19 K; and iii) 
the spinning dust due to t he electric dipole, modelle d as the 
linear combination of two lDraine & Lazarianl d 1998h models 
for warm neutral medium and molecular cloud (0.&5WNM + 
0.35MC). The red line represents the co-added SED model 
of the three contributions. 



Faraday scree n mod el. As mentioned in Sect. [2] 
Reich & Reichl (12009b have recently suggested that 
G159. 6-18.5 is acting as a Faraday screen (FS) hosting 
a strong regular magnetic field, which would rotate the 
background polarized emission, being the responsible 



mecha nism of the signal detected in iBattistelli et al 
(12006b . However, this mechanism produces negligi- 



ble effects on the expected polarized emission at the 
WMAP frequencies. To show this, we can use the 
estimates of the background emission from the last 
two columns in table [3] and take the rotation mea- 
sure RM = 190 rad rrT 2 as a face value to evaluate 
the effect. Being conservative, we can assume that 
the full background emission is rotated by the FS. In 
that case, for the relevant rotation angles (computed as 
9 - RMA 2 ), we obtain that the expected polarization 
signal is smaller than 0.1 Jy (or 0.2%) at 23 GHz, and 
becomes even smaller at higher frequencies (Ka, Q and 
V). Therefore, we can safely neglect its contribution at 
the WMAP frequency bands. 

Implications on the anomalous emission models. 

Figure [6] illustrates the implications of our measure- 
ments on the existing models of ano malous emission 



The fi gure also includes the value of IBattistelli et al. 
d2006h . because even in the case that this result only 
provide an upper limit to the polarization of the 
anomalous emission, it still is a very valuable con- 
straint. 

In Figure [6] we have considered two families of the 
so-called spinning dust models: i) the electric dipole 
emission and the resonance relaxation ( hereafter ED) 
proposed by lLazarian & Draine and ii) the 

magnetic dipole emission (hereafter MD) proposed in 
Draine & Lazarianl (1 1999b . For this second case, the 
number of variants is very large, so here we will re- 
strict ourselves to those models where the grains that 
dominate the 10-100 GHz emission consist of a single 
magnetic domain. These models were computed by 
Draine & Laza rian d 1999b for metallic Fe, and for an 
hypothetical material X4 (defined by them). Both re- 
sults will be used here, for three different grain shapes, 
namely 1 : 1.25 : 1.5 (dotted lines), 1 : 1.5 : 1.5 
(dashed lines) and 1:2:2 (dot-an d-dashed lines). For 
further details on these models, see lDraine & Lazarian 
(Il999b . 

We note that in addition to the anomalous emission, 
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Fig. 5. — : Our measurements on the total fractional linear 
polarization IT (black arrows) in the G159.6-18.5 region. We 
also include the resu lt of IT = 3.4*}-*% at 11 GHz from 
iBattistelli etal] d2006h . For illustration, the solid lines show 
the predicted fractional linear polarization due to thermal 
dust in the region for three values of IX/,,,,, namely 6.0% 
(black line), 3.6% (blue line) and 2.0% (red line), which 
are taken to be co nsistent with range of values allowed by 
iKogut et al.l d2007l) . At the frequencies around the peak of 
the anomalous emission, the constraints are of the order of 
1 per cent, and can be directly translated into constraints on 
the polarization of the anomalous emission. 



when plotting the different models in figure [6] we have 
also included the contribution of the polarized thermal 
dust emission, using the value Yl^usi = 6.0%. In prac- 
tise, this is done as described by this equation 



n„, = 



dust 



la + Idust + I// 



(1) 



where the model for the linear polarization of anoma- 
lous and the thermal dust emission can be represented 
by P a = I a n fl and P dust = ld us ,U dllst respectively. The 
\iust, 'iff and Idusi contributions were obtained in the 
previous section. 

From Figure [6] we see that the measurement at 
1 1 GHz can not completely rule out the model of dust 
grains of X4 with semiaxes 1:1.5:1.5. However, when 
including the upper limits derived here from WMAP 
data, all the considered models based on magnetic 
dipole emission can be ruled out. In contrast, all these 
measurements are in agreement with the predictions of 
the polarization fraction fro m the electric dipole and 
resonance relaxation theory ( Lazarian & Drainell2000l) 
at this frequency range (see Fig. [7]). 

5. Conclusions 

A detailed understanding of the physical mecha- 
nism responsible for the anomalous microwave emis- 
sion is still needed. The polarization properties of 
those regions with anomalous emission, as G159.6- 
18.5 , are very useful tools to disentangle among the 
different models proposed in the literature. 

Here, we used the WMAP7 data to study the in- 
tensity and polarization properties of the emission in 
the G159.6-18.5 region. In intensity, our results con- 
firm the presence of the anomalous microwave emis- 
sion. We present an updated SED model that takes 
into account the contribution of the electric dipole 
emission of very small dust grains rapidly rotating 
dDraine & Lazarianlll998l) . 

Concerning the polarization in the region, we 
present the first constraints on the polarization prop- 
erties of the anomalous microwave emission at high 
frequencies (23-94 GHz). Due to the fact that more 
of the 90 per cent of the emission in G159.6-18.5 is 
diffuse and extended ov er scales larger than 40 ar- 
cmin dTibbs et alj |2010|) . the angular resolution of 
WMAP is not a limitation when constraining the po- 
larization deg ree of the emission in the region. Al- 
though we find no detection of polarization in any 
of the five frequency bands, the derived upper lim- 



11 



I 0.0 



1 .0 



0.1 




n (t. dust) = 6.0% 
n (Spinning Dust) : 
Electric Dipole (ED) 



MD X4 (1:1.25: 
MD Ee (1 :1 .25 
MD X4 (1:1,5: 
MD Ee (1:1.5: 
MD X4 (1 
MD Ee (1 



1.5) 
1.5) 
1.5) 
1.5) 
2:2) 
2:2) 



v (GHz) 

Fig. 6. — : Constraints on several anomalous emission mod- 
els based on our results in the G159.6-18.5 region. Data 
points and upper limits have the same meaning as in Fig. [5] 
The lines correspond to the fractional linear polarization for 
different spinning dust m odels: i) the electric dipole from 
lLazarian & Drainel J2000h (b lack solid line); and ii) the mag- 
netic dipole emission from jDraine & Lazariaril 19991) . that 
predicts a different frequency behaviour (blue and red lines) 
depending on the composition and shape of the grains. For 
completeness, we have added to all models the contribution 
of polarized thermal dust emission with a polarization frac- 
tion of Tldusi = 6.0%. Our results exclude the magnetic 
dipole emission models as the physical process responsi- 
ble of the observed polarization in G159.6-18.5 . In con- 
trast, these are consistent with the expected linear polariza- 
tion from the electric dipole emission and resonance relax- 
ation. 



00 



its allow to exclude a significant number of models 
based on the magnetic dip ole emission of dust grains 
(IDraine & Lazarianll 1999b as the physical process re- 
sponsible of the observed polarization. 

The main conclusion is that the combination of our 
constraints in the range 23-94 GHz, with the the mea- 
sureme nt at 1 1 GHz from the COSMOSOMAS exper- 
iment (B attistelli et al. I l2006l) are consistent with the 
expected linear polarization arising from the electric 
dipole model with resonance relaxation. 

Further observations with higher sensitivity and an- 
gular resolution in this and other regions will be valu- 
able to understand the nature of this anomalous emis- 
sion in polarization. 
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Microwave Background Data Analysis (LAMBDA). 
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Fig. 7. — : Constraints on the grain alignment for both reso- 
nance and Davis-Greenstein relaxation models for grains in 
the cold interstellar medium as a functio n of the frequency. 
Both theoretical curves are taken from lLazarian & Draing 
For the resonance relaxation the saturation effects 
are neglected, which means that the upper curves correspond 
to the maximal values allowed by the para magnetic mech- 
anism . The data point corresponds to the iBattistelli et al.l 
result, while the three upper limits are the results in 
this paper. 
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